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Abstract 

A laser interferometric detector of gravitational waves is studied and a complete solution 
(to first order in the metric perturbation) of the coupled Einstein-Maxwell equations with 
appropriate boundary conditions for the light beams is determined. The phase shift, the 
light deflection and the rotation of the polarization axis induced by gravitational waves are 
computed. The results are compared with previous literature, and are shown to hold also 
for detectors which are large in comparison with the gravitational wavelength. 
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1 Introduction 



The detection of gravitational waves is an important goal of contemporary physics and in 
recent years, research in the field has shifted towards laser interferometric detectors ([0-0] 
and references therein). Large laser interferometers are presently under construction and in 
the planning stages in various countries. Remarkable technological progress has been made 
in the effort to detect the weak gravitational waves that are believed to be emitted from 
extragalactic sources. 

The analysis of a laser interferometric detector is usually carried out in Fermi normal 
coordinates (FNC), which are associated with a freely falling observer carrying the detector 
d, [U. In order to introduce these coordinates, one must restrict oneself to considering detec- 
tors with size a much smaller than the reduced wavelength X/2n of the gravitational waves 
to be detected. This condition is satisfied only marginally in some cases jy], and as larger 
interferometers are being planned, it will be satisfied even less and ultimately, not at all in 
certain frequency bands. Most detectors are conceived to operate at gravitational wave fre- 
quencies Vg < 10 KHz, with maximum sensitivity up to 1 KHz. The Italian-French project 
VIRGO plans the construction of an interferometer with arm size a ~ 3 Km, and maximum 
sensitivity around 100 Hz. The detector in its final stage is expected to operate in the range 
from 10 Hz to 1 KHz. For u g ~ 1 KHz, 2ira/X g ~ 0.06, and in the extreme case v g ~ 10 
KHz, 2ira/X g ~ 0.6. Beam detectors in space for the detection of low frequency gravitational 
waves (v g ~ 10 _1 Hz) have also been proposed (|| [7j; ||] and references therein), although the 
necessary technology is not available at present; for these systems, a ~ 10 6 Km and v g ~ 10 _1 
Hz give 2ira/X g ~ 2. Therefore, it appears convenient to generalize known results to the case 
a > \g/2ir. Moreover, since energy localization and transfer in general relativity are not in as 
nearly a well-established state as they are in the rest of physics, it is particularly important 
that the theoretical predictions for the effects of gravitational waves in their interaction with 
proposed detectors be established with great care. The advantage of laser interferometry as a 
detection mechanism is that it requires no reference to an energy transfer mechanism (or even 
whether such a transfer actually occurs ||) but rather, can be phrased entirely in terms of 
phase shifts. In fact, this carries with it the additional advantage of being expressible in terms 
of scalar quantities upon which all observers must agree. In a recent paper, Lobo || usefully 
focussed upon these aspects. 

The formalism for the interference between electromagnetic and gravitational waves was 
established many years ago [1C] and the perturbations to the electromagnetic waves in the 
case of weak fields were solved for longitudinal and transverse orientations. Moreover, they 
were solved with reference to well-defined boundary conditions. In Lobo's work, the phase 
shift is calculated for a wave propagating in one direction and this is added to the phase 
shift for a wave travelling in the opposite direction. However, to correspond to the actual 
physical configuration of a laser interferometer, it is necessary to be more specific and analyze 
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the response to a gravity wave of electromagnetic waves which travel equal distances along the 
interferometer arms at right angles to each other and then reflect at the mirrors which terminate 
these arms. Only when this is done can one have confidence in the expected interference which 
one seeks as the waves are recombined. In this paper, we use the electromagnetically gauge- 
invariant perturbation formalism of Ref. [[H]] in conjunction with the boundary conditions 
appropriate to a laser interferometer, to derive the expected phase shift as the reflected waves 
recombine with interference. Using this formalism, we are also able to compute the extent of 
the deflection of the laser beams due to the interference with the gravitational waves. While 
this is a first order effect in the metric perturbation, it only induces a second order contribution 
to the phase shift. The value of the deflection angle is then compared with the results from 
previous analyses made in a different context. In addition, we show that gravitational waves 
do not induce any rotation in the polarization plane of electromagnetic radiation. 

The plan of the paper is as follows: in Sec. 2 an idealized model of an interferometric 
detector is described, and the Maxwell equations in the presence of a weak, monochromatic, 
optimally polarized gravitational wave are solved for an example case. In Sec. 3, the solu- 
tions for the actual interferometer are derived from this case by means of suitable coordinate 
transformations, and the phase shift in the light coming from the two interferometer arms is 
computed. In Sec. 4, we take into account the deflection of light beams and changes in the 
polarization induced by gravitational waves, showing that these effects are negligible. Finally, 
the phase shift is compared with previous literature in Sec. 5. 

2 An idealized beam detector 

We consider a simplified model of a laser interferometric detector, whose geometry is given 
in Fig. 1. The detector consists of a Michelson interferometer, in which three masses are 
hanging in the plane z = 0. The periods of the suspension pendula are much greater than the 
gravitational wave period T g , so the masses can be regarded as freely falling. A beam splitter 
at the origin splits the light from a laser lamp into two beams of equal intensity which travel 
along the x and y axes, and are reflected by mirrors placed at a distance a from the origin, 
and perpendicular to the x and y axes, respectively. The reflected beams are collected at a 
photodiode which compares them and detects any differential effect induced by an impinging 
gravitational wave in the interferometer arms. 

We assume that (see Ref. || for a discussion of these assumptions) 

• the electromagnetic field in the laser beam can be treated as a test field; 

• the mirrors at x = a, y = a are perfectly reflecting; 

• the transverse dimension of the light beam is much smaller than the gravitational wave- 
length X g ; 
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the light frequency is much larger than the gravitational wave frequency, 

LO » Ug . 



(2.1) 



We consider a monochromatic gravitational wave travelling along the positive z axis with 

a single polarization, and with polarization axes along the "L" of the interferometer. The 
spacetime metric is given byQ 

Qiw = V + hfj, u (2.2) 

in an inertial frame of the Minkowskian background metric rjuu, and \huu\ << 1. For a polarized 
beam in the TT gauge, the only nonvanishing huv describing the gravitational wave are 
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where A is a (small) constant and k g = uj 9 /c. Although we will restrict our attention to 
this particular direction of propagation and to this polarization of the gravitational wave, the 
calculation can be generalized to arbitrary incidence and polarization, in the manner of Ref. ||. 

The unperturbed electromagnetic field in the x arm of the interferometer has the only 
nonvanishing components 
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where Eq =constant, k = u/c, and where only the real components of the complex fields need 
to be retained. The unperturbed field in the y arm is given by 
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These fields satisfy the flat space Maxwell equations 

F(°^% = 0, (2.8) 

, (2.9) 

with perfectly reflecting boundary conditions at the mirrors located at (a, 0, 0) and (0, a, 0) 
(see Fig. 1). We decompose the electromagnetic field tensor as the sum of the unperturbed 
field and a small perturbation induced by the gravitational wave 



p(o) , p(0) 

fj.u,p 1 M Vp,fl 



F (0) 



(J.V 



(2.10) 



x The metric signature is +2, and c denotes the velocity of light in vacuum. Greek indices run from to 3 
and Latin indices run from 1 to 3. We will perform computations to first order in the metric perturbation h^ v . 
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with \F$)\ << \Ffu/\. The field equations for the perturbations F^' in the TT gauge have 
been obtained in Ref. [1C]: 
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(2.11) 
(2.12) 



It has been shown in Ref. [11] that other expressions in the literature are incorrect because 
they do not account properly for the Lorentz condition^. Equations ( |2.11| ) and ( 2.12j ) must be 
solved in both arms using perfectly reflecting boundary conditions. We solve these equations 
in a special orientation which does not correspond to an actual interferometer arm in our 
geometry depicted in fig. 1. We then obtain the solution to Eqs. ( 2.11] ) and ( 2.12j ) for the 
actual interferometer arms by means of suitable coordinate transformations. The fictitious 
system is composed of an electromagnetic wave propagating along the z axis, from z = to 
z = a, where it is reflected by a perfect mirror; the system is perturbed by a gravitational wave 
described by 

/in = -h 33 = Acos(k g y - u g t) . (2.13) 
The unperturbed electromagnetic field is given by 
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Equations ( [2.11 ) for the electromagnetic tensor perturbations are 
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(2.15) 

(2.16) 
(2.17) 
(2.18) 
(2.19) 



from which we can easily derive an inhomogeneous wave equation for F^ 
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) and (^.20|) 
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are linear both in h^ v and in , but not in these two kinds of 
variables simultaneously. As a consequence, one cannot use the complex representations for both fields at the 
same time, retaining only their real parts in the final results. However, one is always allowed to make use of the 
identity cos 2 = (e iz + e~") /2. 
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where 
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The solution of Eqs. ( |2.16 )-( 2.19| ) and ( 2.1 2| ) with the perfectly reflecting boundary condition 
at z = a and to first order in h is 
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and the other F^) = 0. We can now use this example case to write the solutions to Eqs. (|2.11 ) 
and (2.12) for the actual interferometer arms. 



3 Solutions of the Maxwell equations and computation of the 
phase shift 



x arm: 

By means of the coordinate transformation {x^} i— > {x'^}, with 
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the unperturbed electromagnetic field tensor given by Eqs. fl2.4|) and ( |2.5| ) transforms in such 
a way that its only nonvanishing components in the {x ,f1 } system are 
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The metric perturbation becomes 

h'n = -h'33 
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Equations Q, (|J) and (|J) coincide with Eqs. (|T4|), ( $1$ ) and (|2l3|) , the solutions of 
which are given by Eqs. ( |2.23[ )-( p.25 ). Transforming back to the {x^"} system, these become 
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and the other F^J = 0. 
y arm: 

We now perform the coordinate transformation 
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Thus, we reproduce the case described by Eqs. ( |2.13| )-( |2~T5| ), provided we make the following 
substitutions: 
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The electromagnetic field perturbations are then given by Eqs. (g2g)-(g2g), with these substi- 
tutions. Transforming back to the {x^} coordinate system, we find that the only nonvanishing 
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In the TT gauge, the effect of the gravitational wave is to generate two sideband components 
of the electromagnetic signal at frequencies Q^, which propagate together with the carrier 
of frequency to. This can be seen as a phase shift at the photodiode, where the two signals 
reflected by the mirrors are compared. 

In the x arm, in the limit to >> u> g , the backward propagating component of the electric 
field on the z = plane is 
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and comparing Eqs. (3.17) and ( |3,18| ) at the arrival time r of the signals at the photodiode (r 
is the round trip time for the laser light in each of the interferometer arms), we find the phase 
shift for the x arm: 

5(j> x = — — sm(u g T) . (3.19) 

I UJg 

Analogously, we obtain the phase shift for the y arm 



A LO 



- — sm[fj g T) 



(3.20) 



which is equal and opposite to 6<j> x . The differential phase shift between the signals coming 
from the two interferometer arms is thus 
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expressed in a form which is convenient for the purposes of comparison. 



4 Beam deflection and polarization shift 

We now consider the deflection of the laser beams induced by their interaction with the grav- 
itational waves. The fact that gravitational waves curve light rays propagating through them 
is well known, and has been studied in the geometric optics approximation by various authors, 
mainly in view of its possible astrophysical consequences [|l|-[|ll|. To the best of our knowl- 
edge, it has not been discussed previously in the context of interferometry. In principle, it 
could be of significance because it implies a change in the distance of beam traversal and as a 
consequence, a variation in phase of the waves. However, we will demonstrate that this phase 
shift is negligible in comparison with that discussed in Sec. 3 under the normal conditions 
currently considered. 

We now analyze the configuration described by Eqs. ( |2.13| )- (l2.15| ): The unperturbed Poynt- 
ing vector of the electromagnetic waves becomes 

5(0) = c £(o) x £(o) = c E (S» H (o) g z (4 J) 

(where et denotes the unit vector parallel to the i-th axis). In the presence of the gravitational 
wave, the Poynting vector is 

g{tot) = CexH = + S« , (4.2) 

where 

= h [(^W + E * )H f ] ) ^ ~ ^^H] + 0(2). (4.3) 
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Clearly, the vector S^ tot ' is rotated in the (y,z) plane with respect to S^°\ We consider, for 
the sake of simplicity, only the forward propagating component of the electromagnetic wave. 
Let 9 be the angle between the vectors S^ rward and S^ ward in the (y, z) plane where they 
lie, i.e. the deflection angle. We have 



n{tot) 
y forward 



forward 



sm/ 



from which we obtain (we omit the subscript "forward" in the following) 



(1) 



sm( 



\H, 



(0) 
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In the limit to » uj g , in the plane y = 0, and expanding sin 9 to first order, we find 

9 = — cos(uj g t) . 



(4.6) 

Note that 9 = 0(h). Equation ( [4.6[ ) can be compared with the deflection angle obtained in 
previous analyses performed in the geometric optics approximation. If we consider a photon 
whose unperturbed path is parallel to the z axis, and with four-momentum 
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where 5p^ are small deflections, the equation of null geodesies give 
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(i = 1,2), where the integral is computed along the unperturbed photon path from the light 
source S to the observer O. This formula can be found, e.g., in Refs. |17]-[20|. In our case we 
obtain, from Eq. Q2.13 ), 



5p x 
8p y 



0{2) 

4cos(<V)+0(2) , 



(4.9) 
(4.10) 



which confirm that the deflection takes place in the (y, z) plane, and the value of 9 is as in 
Eq. (4.6). This provides a check of our results from an independent source. It is of interest to 
consider the additional phase shift at the photodiode due to the change in path length travelled 
by the photons under the deflection. If 5li is the extra length for the i-th arm (i = 1,2), the 
phase shift it induces is 5<pi = 2ir5li/\. It is immediate to see that, if li is the total length 
travelled by the photon, 



SL = L 



COS I 



(4.11) 
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The additional phase shift introduced by the deflection of light is thus Sfa = 0(h 2 ), and is 
completely negligible. This is true also for the differential phase shift in the two arms of the 
interferometer. 

An interesting property of the interaction between gravitational and electromagnetic waves 
can be seen from the solutions ( |2,23| )-( |2,25| ) of the coupled Einstein-Maxwell equations. We 
note that the unperturbed electric and magnetic fields are polarized along the x and y axes, 
respectively. The perturbation to the electric field is parallel to the x axis, while that of 
the magnetic field has components along both the y and z axis. This is consistent with the 
result [21 1 that a weak gravitational wave induces no rotation in the plane of polarization of 
electromagnetic radiation, to first order in h. In fact, let us consider a given point O on the 
light beam, at a given time, and the two planes V\ and V2 passing through O, which are 
defined by the pairs of vectors (E^,H^) and {E^ tot \ H^), respectively (see Fig. 2). The 
normals to these planes are parallel to the Poynting vectors and S( tat > respectively. V\ is 
associated with the forward propagating electromagnetic wave at O at the given time, in the 
unperturbed case (no gravitational waves) and V2 with that corresponding to the perturbed 
wave. V\ and V2 intersect along their common x axis. The plane V2 can be obtained by 
rotating V\ by an angle 6 around the x axis. The projection of the magnetic field e V2 

on V\ lies on the direction of (i.e. on the y axis of Pi). In this sense, the variation of the 
direction of the magnetic field due to the gravitational wave arises solely from the deflection of 
the light beam. No rotation of the electric and magnetic field, with respect to the unperturbed 
case, takes place in a plane orthogonal to the Poynting vector (in the sense of the background 
Minkowski metric). Although these considerations could be made more rigorous by using 
advanced concepts from differential geometry, it is unnecessary, since a general proof of the 
fact that gravitational waves do not rotate the polarization vector of electromagnetic waves 
propagating through them, to first order, can be found in Ref. [21]. 



5 Discussion and conclusion 

In Sec. 3 we have presented a complete solution (to first order in h) to the Maxwell equations in 
the field of a weak gravitational wave, taking into account the boundary conditions. Our result 
in Eq. ( |3.21| ) disagrees with that by Lobo by the presence of the factor cos(o; 9 r/2). This factor 
is absent in Lobo's paper ||, but it does appear in a paper by Meers (T^[. This is contrary 
to Lobo's claimed agreement with Meers. To reproduce Lobo's result in our calculations, it 
would require a time interval (0,r/2) to be considered twice instead of correctly considering 
the two distinct time intervals (0, r/2) and (r/2,r). The former approach would give a factor 
2 sin(d; 9 r/2) in our Eqs. ( |3.19D and (|3.20| ), and corresponds to computing the phase shift for 
the forward propagating signal during the time interval (0, r/2), and incorrectly doubling the 
result in order to obtain the phase shift for the reflected signal travelling from the mirror to 
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the photodiode during the time interval (r/2,r). Apparently, the final formulas for the phase 
shift in Ref. B are obtained in this way. The correct computation gives a factor sin(a; 9 r) = 
2 sin(u; 9 r/2) cos(a; 9 r/2) instead of 2sin(u; 5 r/2) in Eqs. ( |3.19| ) and ( p. 20 ). However, it should 
be noted that the distinction is irrelevant in the limit r/T g << 1, wherein both Eq. ( 3.21| ) and 
the result of Ref. || agree with Ref. Q. This limit is not unphysical: if r/T g ~ 1, Eq. ( p. 21 ) 
gives 5(j) ~ 0, and the detector is rendered useless due to the fact that the metric perturbation 
reverses its sign and the phase shift accumulated in the first half period of the gravitational 
wave is destroyed during the second half [|l|. However, the inherent broadband nature of the 
detector makes it sensitive to gravitational wave pulses which are not monochromatic, but 
have Fourier components in a band of frequencies Af 9 , some of which may have co g T less than 



1, but not very small. In such cases the full Eq. (3.21) applies. 



Perhaps, the most promising sources of gravitational waves detectable in the near future 
with beam detectors are coalescing binaries containing compact objects. These systems give 
rise to a characteristic signal ( "chirp" ) with frequency increasing up to 1 KHz in the final stages 



of coalescence. The final moments of the process are not completely understood |U|, 14, 15], 
and might well give rise to higher frequencies not satisfying uj g r << 1. 

In Sec. 4, we used the formalism in Ref. ]l0| to deduce the deflection of the laser beams in 
their interaction with the gravitational waves. We found that, while this is a first order effect 
in h, it only induces a second order contribution to the phase shift. 

Finally, we note that if perfectly reflecting mirrors are placed at x = 0, y = normal 
to the x and y axes, the wave number of the unperturbed electromagnetic field in the cavity 
so constructed is restricted to the values k = nir/a, where n is an integer. This system can 
provide an idealized description of a laser interferometric detector operating with Fabry-Perot 
cavities in its arms and references therein). 
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Figure captions: 

Figure 1: The gravitational wave propagates along the positive z axis; electromagnetic waves 
starting from the beam splitter at the origin are reflected by mirrors placed at x = a, y = a. 

Figure 2: The effects induced by the gravitational wave on the plane of polarization of a 
forward propagating electromagnetic wave. EJ(°\ and are the electric and magnetic 
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fields, and the Poynting vector, respectively, in the unperturbed case. & tot \ #(* *) and 
are the corresponding quantities in the presence of gravitational waves. is the deflection 
angle of the beam. The planes V\ and V2 are the polarization planes of the electromagnetic 
radiation in the unperturbed and perturbed case. 



13 



